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Abstract. The Wolf-Rayet binary system 7^ 
the closest known Wolf-Rayet (WR) star. Rec 
tance has been redetermined by parallax 
the HIPPARCOS astrometric satellite yielding 
nificantly lower than previous estimates (300 - 
Rayet stars have been proposed as a major source 
tic ^''Al observed at 1.809 MeV. The gamma 
COMPTEL has previously reported 1.8 MeV 
the Vela region, yet located closer to the Galac 
the position of 7^ Vel. We derive an upper 1.8 
of 1.1 10-5 7cm-2 s-i (2 a) for the WR star 
distance estimate, COMPTEL measurements pi 
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1. Introduction 

The 1 .809 MeV gamma-ray line from radioactive 
[mean Hfetime (1.07 ± 0.04) 10^ yr (|Endt 19' 



cent nucleosynthesis in the Galaxy. This was the 
ray line detected from the interstellar medium 
1984), and had been predicted from nucleosynt^i 
tions of explosive carbon burning in core-collap 
( |Ramaty & Lingenfelter 1977| ; [Arnett 1977| ). Olher 
duction sites have been proposed as well, 
range of densities (0.5 - 3 10^ g cm~^), temperatures 
3 10^ K), and time scales (1 - 10^^ s) at which j 



g a wide 
(3 10^- 
roton capture 
-''Al: explo- 
sive hydrogen burning in novae and supernovae (jArnould et al. 



on Mg (within the Mg-Al chain) would create 



19800 

(e.g. 

thesi 

drogi 

andl 

of th 



Woosley & Weaver 1995), neutrino-induced nucleosyn- 



i in supernovae ( Woosley et al. 199C ), convective core hy - 
burning in very massive stars (Dearborn & Blake 1985), 

ydrogen shell burning or "hot bottom burning" at the base 
convective stellar envelope in asymptotic giant branch 
(AGlB) stars ( |N0rgaard 1980| ; [Forestini et al. 1991| ). ^''Al nu- 
cleos ynthesis and observations have been reviewed by Clayton 

isi ng (1987);|MacPherson et al. (1995] ); Prantzos & Diehl 
(199^); piehl&Timmes(1998D . 

urrent theories still predict significant amounts of ^^Al 

core-collapse supernovae (Woosley & Weaver 1995 ; 



C 

from 



Thie: 



emann et al. 199q; Timmes et al. 1995 ), the wind phases 



of 



esis - 
: individual 



ties 
mas 



for a 

2«A1 



e supernovae 
^''Al pro- 



low- 
26 Al 
prod 



neon burning in the presupernova and supernova stage 



velry massive stars where ^^Al is produced on the main se- 
quen;e an d ejected mainly in the Wolf-Rayet stage (Langer 
et al. 1995; Meynetetal. 1997 ) or possibly (in case of fast rotat - 
ing s ars) in a red supergiant stage as w ell (Langer et al. 1997 ), 
and irom the most massive AGB stars ( Bazan et al. 1993 ). The 
expe ;ted ^''Al contribution of chemically enri ched novae low- 
ered ifter major revisions of key reaction rates ( Foseet al. 1997 ; 



in 



■^ 



[Starr leld et al. 1998 ). All models suffer from large uncertain 



Al yields, ranging from factors of three to orders of 

ilitude for the various proposed astrophysical sites. 

F ecent results from the COMPTEL telescope aboard the 

Com 3ton Gamma-Ray Observatory showed growing evidence 

young (massive) population dominating the galaxy-wide 

production dDiehl et al. 1995b 



decay of '^^ Al 
0)] traces re- 
first gamma- ^ ^__^^^^.^^^^ 

(^ lahoney et al . jOber lacketal. 1996| ; pbei-lack 1997 
esis calcula 



Knodlseder et al. 1996b; 



Kn6dlsederl997). While 



mass AGB stars and novae can be ruled out as the main 
source, a distinction between supernovae and hydrostatic 
iction in very massive stars appears difficult due to the 
similir evolutionary time scales involved.]] Therefore, detec- 
tion )f individual objects would be essential as calibrator, but 
the s ;nsitivity of current instruments restricts this approach to 
very few objects. Upper limits for five individual super nova 
remnants have been derived by Knodlseder et al. (1996a ) and 
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See Knodlseder (1995) for arguments favouring WR stars as the 
dominant source of Galactic ^*A1. 
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the possibility of interpreting 1.8 MeV emission from the Vela 
region with individual objects has been discussed by Oberlack 



et al. (1994) and Diehl et al. (1999) 



7^ Vel (WR 11, van der Hucht et al. 1981) is a WC8+08 



8.5IIIbinary system at (/, b) = (262.8°, -7.7°), containing the 
nearest Wolf-Rayet star to the Sun at a distance of 258^3]^ pc, 
as determined by parallax measurements with the HIPPAR- 
COS satelhte ( |van der Hucht et al. 1997^ [Schaerer et al. 1997[ ). 
[Another recent investigation determines a spectral type of 
07.5 for the O star ( pe Marco & Schmutz 1999| ).] A previous 



1.8 MeV flux limit(2 cr) for 7^ Vel of 1.9 10'^ 7 cm"^ s"^ had 



been determined by piehl et al. (1995a ), based on observations 
of the first 2 1/2 CGRO mission years. 

In Sect. 2, we describe our search for 1.8 MeV emission 
from 7^ Vel and the derivation of upper flux limits for several 
emission models. Sect. 3 discusses the initial mass range of the 
WR star, implications of our flux limit for stellar models, and 
proposes an alternative interpretation of the "IRAS Vela shell". 
We summarize in Sect. 4. 

2. Data analysis and source models 

The COMPTEL telescope spans an energy range from 0.75 to 
30 MeV with spectral resolution of 8% (FWHM) at 1.8 MeV, 
and performs imaging with an angular resolution of 3.8° 
FWHM at 1.8 MeV within a '^ 1 sr field of view. It features 
the Compton scattering detection principle through a coinci- 
dence measurement in two detector planes (see Schonfelder 
et al. (1993) for details). Imaging analysis and model fitting 
occurs in a three-dimensional dataspace consisting of angles 
(x, fp, <f) describing the scattered photon's direction and the es- 
timated Compton scatter angle, respectively ( piehl etal. 1992 ). 

In this paper, we concentrate on fitting sky models, con- 
volved with the instrumental response, in the imaging datas- 
pace. The present analysis makes use of all data from obser- 
vations 0.1 - 522.5 which have been combined in a full-sky 
dataset, comprising 5 years of observing time between May 
1991 and June 1996. Events were collected within a 200 keV 
wide energy window from (1.7 - 1.9) MeV into the imaging 
dataspace with 1° x 1° binning in (x, ^) (in galactic coordi- 
nates) and 2° binning in (p. The dataspace has been restricted 
to I G [185°, 320°] and |6| < 50° to concentrate on emission 
from the Vela region, but not to loose information from the 
up to 100° wide response cone (at (p = 50°). The instrumen- 
tal and celestial continuum background have been modelled by 
interpolation from adjacent energies, with enhancements from 
Monte Carlo modelling of identified activation background 
components. A previous version of the background handling 
and event s elections have been desc ribed by Oberlack et al. 
(1996) and |<^nodlseder et al. (1996c ), more details on the re- 
cently improve d background ha ndling and the complete dataset 
are reported in Oberlack (1997 ) [Oberlack et al., in prep.]. 

For the derivation of upper limits, the maximum likelihood 
ratio test has been applied ( [Cash 1979 ): A null hypothesis Hq 
is compared with an extended alternative hypothesis Hi, which 
includes q additional continuous model parameters, using the 



likelihood function C, which is the product of the Poisson prob- 
abilities pk in N dataspace cells: 



N 

c = XIp. 

fe=l 



Pk 



^ e-'^" for Aifc > 



for ^fe ^ 0, nfc == 
for ^fe == 0, rife > 



(1) 



where n^. is the number of counts in cell k and 



Mfc 



s=l 



I, bgd 



(2) 



is the predicted number of counts due to sources and the back- 
ground model, which includes the scaling parameters a^'^\ b 
varied in the fit. Each source model (s) can be described by 
a (normalized) flux map {/„■ } convolved with the response 
matrix Rjk'- 



(s) 



:rm 



(s) 



(3) 



If the null hypothesis is true, the probability distribution of 
the ratio of the maximum likelihood £1 achieved by fitting Hi 
to the data over the maximum likelihood Cq achieved by fitting 
Hq to the data can be described analytically: 



p 2 In 




-pixD 



(4) 



where p(Xg) is the tabulated x^ probabiUty distribution with q 
degrees of freedom. 

Fig. |l] shows a Maximum-Entropy deconvolved map of the 
1.8 MeV emission from the Carina / Vela / Puppis regions. 
While no 1.8 MeV flux is detected from the position of 7^ Vel, 
significant extended emission is observed in nearby regions, 
with an intensity peak around (267?5, — 0?5). Due to the broad 
COMPTEL response, such emission needs to be modelled. We 
test its impact on flux limits for 7^ Vel with four different emis- 
sion models (additional to the background model), guided by 
the deconvolved image, by candidate ^^Al sources in the region, 
and by results on the galaxy-wide ^^Al distribution. These are 
the tested source models in addition to the background model: 

a. A single point source at the position of •y'^ Vel. 

b. A detailed model describing the observed emission em- 
pirically and including 1.8 MeV candidate sources of the 
region (Fig. nb: an exponential disk with emissivity oc 
cxp{— i?/i?o} • exp{ — [zj/zo}, galactocentric scale length 
i?o = 4 kpc, and scale height zo = 180 pc as an approx- 
imation to the large-scale ^^Al distribution, an additional 
homogenous stripe around I ~ 310° for simplified model- 
ing of excess emission in this region and 8 point sources 
including 7^ Vel. Details are listed in Table IT]. 

c. An intermediate model: 2 point sources at the positions 
of maximum intensity in Fig. |I] at (/,6): (286?0,0?0), 
(267?5, -0?5) plus a point source for 7^ Vel. 
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Fig. 1. COMPTEL 1.8 MeV map of the Carina / Vela / Puppis 
regions based on a deconvolution using the Maximum Entropy 
Method (pberlack 1997). Contour lines are equally spaced 



with steps of 2.9 10^^ 7cm^^ s~^ sr~^. The dashed line in 
dicates the IRAS Vela shell, the marked point source models 
are listed in Table [l|. 



Exponential disk: 



Ro=4 kpc, zo = 180 



Homogenous stripe: 300° < I < 320° \b\ < 3° 



Point sources {l,b) 
*i262°S,-7°7) j^Yel + (238?0, 4?0) Puppis feature 

D (256?0,-4?7) C Pup o (263?6 -2?8) Vela pulsar 

A (266?2 -1?2) SNR "Vela Junior"^ + (267?5,-0?5)Vela feature 
A (272?2,-3?2) SNR G272.2 - 3.3 + (286?0, 0?0) Carina feature 



Table 1. Source components of the detailed model (b). For the 
exponential disk, scale length and scale height are given as well 
as the range of the homogenous stripe. Point sources reflect 
candidate ^^Al sources in the reg ion and additional emission 
features appearing in the map. ^ Aschenbach (199q ); lyudin 
etaL(1998) 

d. Like model (c) but replacing the point source for 7^ Vel 
by a model for the IRAS Vela Shell: a spherical shell 
with a thickn ess of 10% of the 8° radius around {l,b) = 
(263°, -7°) ( [Sahu & Blaauw 1993| ). 



The first three models assume that all ^^Al from 7^ Vel is 
kept within an observed 'ejecta-type' wind shell around th e 
binary system with a diameter of 57' (Marston et al. 1994). 
Given COMPTEL's angular resolution, this can be modelled 
by a point source. For completeness, we consider a fourth 
model, which places ^^Al into the so-called IRAS Vela shell, 
an extended (^ 8° radius) structure almost centered on j'^ Vel, 
which has been identified by [Sahu (1992 ) from IRAS 25/^m 
- 100/im maps and investigations of cometary globules. This 
is a region of relatively bright Ha emission within the Gum 



nebula (Chanot & Sivan 1983), where young stellar objects are 
forming ( ^rusti et al. 1992 ). Sahu found the Vela shell to be a 
structure distinct from the Gum nebula and presented an inter- 
pretation as supershell from the aged Vela OB2 association of 
which 7^ Vel may be a member ( de Zeeuw et al. 1997| ). (As a 
consequence of HIPPARCOS parallaxes, 7" Vel would be lo- 
cated on the very near side of the association.) In this case, part 
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Fig. 2. Determination of upper flux limits with the maximum 
likelihood ratio test using 4 different emission models (see 
text). 

of the ejected -^Al may have traversed the gas bubble and be 
accumulated in the dense outer shell. We discuss a different 



interpretation, however, in Section 3.1 . 

Fig. shows the logarithmic maximum likelihood ratio 
2 \n{C/Ca) as a function of model flux. No significant flux 
from 7^ Vel is found for any of the tested models. The follow- 
ing 2 (T upper limits are derived: 



Model a: 


/ 


< 


1.6 10- 


"^ 7cm^^s~^ 


Model b: 


/ 


< 


0.9 10- 


"^ 7cm^^s^-^ 


Model c: 


/ 


< 


1.1 10" 


'^ 7cm^^s^^ 


Model d: 


/ 


< 


1.9 10" 


"^ 7cm~^s"^ 



Given our 1.8 MeV map of the region, model (a) seems overly 
simplistic. Lacking alternative source models, part of the ob- 
served Vela emission is attributed to 7^ Vel and yields the most 
conservative upper limit for a point source model. Model (b) 
contains many free parameters and may approach an over-fit of 
the data, where statistical fluctuations of the background are fit- 
ted. Although the fitted background scaling factor is lowest for 
this model, the lowest (even slightly negative) flux is attributed 
to 7^ Vel. Consistent with the map in Fig. [I], the flux from Vela- 
Puppis is better described by sources within the galactic plane 
than by 7^ Vel. Model (c) is "intermediate" in that it accounts 
for the strongest features in the map with a minimum of free 
model parameters. We adopt its result for the 7^ Vel flux limit 
and consider models (a) and (b) as the extreme values for the 
systematic uncertainty due to the choice in modelling of other 
emission close (~ 5° - 10°) to 7^ Vel. The limit from model (d) 
is highest due to the large extent of this source model. Yet, we 
do not consider this model a likely representation of ^^ Al from 
7^ Vel, even if the structure itself may well be related to the 
binary system, as discussed in the next section. 

The 1.8 MeV flux direcdy translates into the "ahve" ^^Al 
mass in the circumstellar medium, for a point source via: 



/i. 



8 MeV 



= 1.8 10"''7cm"^s- 



258 \ ^ M26 [Mg 



rf[pc] 



1.0 10- 
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Therefore, our 2 a upper flux limit corresponds to a maximum 
^^Al mass from 7^ Vel of 



Af^^Rii < (6.3t?i) 10-5 M 







(5) 



where the la distance uncertainties from the HIPPARCOS 
measurement have been taken into account. 



3. Discussion 

3.1. Interpretation of the IRAS Vela shell 

Is our upper limit for the ^^Al yield from 7^ Vel realistic, being 
derived for a point source model, or should we rather consider 
the higher value implied by potential ^^Al accumulation in the 
IRAS Vela shell? 

The refined distance of 7^ Vel suggests a new interpreta- 
tion of the IRAS Vela shell as the main sequence bubble of the 
WR progenitor star in the 7^ Vel system. This would make a 
significant ^^Al contamination of the shell unlikely since this 
isotope is expected to appear at the stellar surface only in later 
stages of stellar evolution, together with the products of core 
hydrogen burning, after the hydrogen shell has been expelled. 
Evidence for our interpretation stems from observations of in- 
terstellar reddening by [Franco (1990 ), who studied two regions 
within the Gum nebula, one of which in projection to the IRAS 
shell. Only this field showed clear evidence for a dust wall at a 
distance of 200 ± 20 pc, interpreted by the authors as the near 
edge of the Gum nebula. If this is now interpreted as the near 
edge of the IRAS shell, its angular extent would correspond to 
a radius of 32 pc and the centre would be placed at a distance 
of 230 ± 20 pc, well within the la uncertainty of the 7^ Vel 
distance. This would argue against a supershell interpretation 
since the distance to the centre of Vela 0B2, the assumed origin 
of the supershell, has been precisely measured by HIPPARCOS 
to 415± 10 pc (de Zeeuw et al. 1997). Scaling down the (uncer- 
tain) mass estimate for the shell by Sahu from the distance of 
450 pc she assumed to the reduced distance of 7^ Vel, yields a 
mass of ~ 2 X 10^ M©. Combined with the observed expansion 
velocity of ^ 10 km s^^ (Sahu & Sahu 1993) this corresponds 
to a kinetic energy of ^ 2 x 10'''^ J (2 x 10^° erg), within the 
range of stellar wind energy release by a massive star. In a hy- 
drodynamic model coupled with a stellar evolution code for a 
60 M0 star, 3arcia-Segura et al. (1996) even find a total en- 



TIT 



ergy release of 3.3 x lO"*"* J (70% - 80% thereof ejected before 
a "luminous blue variable" [LBV] phase) with a 45 pc radius 
O-star bubble. This number only scales weakly with ambient 



density (ex n' 



-1/5 



20 cm ■^ assumed in the model), but 



is considered an upper limit by the authors because effects like 
heat conduction or cloud evaporation are ignored here. Within 
the uncertainties of the model, our interpretation of the IRAS 
Vela shell seems therefore plausible. For the further discussion 
we will hence adopt our 1.8 MeV mass limit which was de- 
rived from the assumption that all ^* Al from 7^ Vel is contained 
within a region of ~ 1° in diameter around 7^ Vel. 



3.2. Current mass estimate for 7^ Vel 

Model predictions of ^''Al yields for massive stars are a strong 
function of initial stellar mass, e.g., Af26 c< M?-^ (Meynet et al. 
1997). A direct determination of the WR star initial mass (and 
thus its predicted ^''Al yield) from comparison of luminosity 
and effective temperature with stellar evolution tracks, as is 
typically done for other types of stars, is not feasible due to the 
lack of sufficiently accurate models of WR wind atmospheres 
(and the additional complications from the binary nature of this 
stellar system). While binarity makes modeling more compli- 
cated due to additional degrees of freedom in parameter space it 
allows measurement of current masses, which can be matched 
with theoretical predictions together with the generic spectral 
type of the WR star Spectroscopic determination of Mi. 2 sin'^ i 
(where i is the inclination) based on Doppler-shifted absorp- 
tion (O star) and emission lines (WR star) led to contradictory 
results (fPike et al. 19831; |Moffat et al. 19861). A recent redeter- 



mination of orbital parameters by Schmutz et al. (1997) yields 
spectroscopic masses of: 



MwR sin^ i = 6.8 ± 0.6 Mq 



Mosm^i = 21.6 ±1.1 



Mr. 



They reject an earlier inclination measurement from pola risa- 
tion data of i = 10° ± 10° by |St.-Louis et al. (1988| ) and 
rather state wider inclination limits of 57° < i < 86.3° from 
other evidence, corresponding to a factor l/sin"^i = 1.0 - 
1.7 or a range from 6 to 12 M© for the current mass of the 
WR star Relying on a mass-lumin osity relation for the O star 
from single star evolution models, [Schaerer et al. (1997 ) de- 
rive Mo = 29 ± 4 Mq, which, in turn, yields a consistent, 
but model-dependent, inclination estimate of i = 65° ± 8° 
or a mass estimate f or the WR star of Mwr = Sl^^ Mq 
( Schmutz et al. 1997 ). A different analysis of the same spectral 
data yields a slightly higher luminosity for the O star, leading 
to a slightly higher O star mass estimate of 30 ± 2 M© (De 
Marco & Schmutz 1999), yet, with the same mass-luminosity 
relation from single-star models. 

Another observational hint on the total mass of the binary 
system has been derived from interferometric measurements of 
the major half-axis of the binary system a" = (4.3 ± 0.5) mas 
by Hanbury Brown et al. (1970). Kepler's law leads to a con- 
sistent, but barely constraining, mass estimate: 



MwR + Mo = 



(2^)^ (g" ■ df 

rp2 



G 



(30 ± 16) Me 



(6) 



The large relative uncertainty in total mass of ~ 54% is due 
to the uncertainties in a" and d to about equal amounts. Addi- 
tional systematic errors, however, may affect the determination 
of a" by a fit in which other quantities like orbital parame- 
ters, inclination, and brightness ratio of the two stars had to be 
taken as fixed parameters due to limited data quality. Those val- 
ues were taken from other measurements available at the time. 
Most notably, an eccentricity of 0.17 had been assumed which 
is distinctl y lower than a recent value of 0.326 ±0.01 (S chmutz 
et al. 1997). Future interferometric measurements could im- 
prove this situation significantly. 
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3.3. Single star models 

Given the remaining uncertainties in current mass determina- 
tions, not surprisingly initial mass estimates and predicted ^^Al 
yields vary significantly. To start with the simpler case, we first 
discuss our results in the context of single star models, assum- 
ing that the stellar structure of the primary (defined in terms of 
which star evolved first, i.e. now the WR star) has not been sig- 
nificantly altered due to the presence of its binary companion 
and that all of the mass expelled from the primary has reached 
the ISM without being captured by the secondary. This corre- 
sponds to fully non-conservative mass transfer, which is usu- 
ally parametrized by the fraction /3 of mass that is accreted by 
the secondary out of the mass lost by the primary (i.e. fully 
non-conservative mass transfer means /? = 0). For the Geneva 
single star models, a minimum initial mass of about 40 M© is 
requ ired at solar metall icities to yield a WC-type Wolf-Rayet 
star( Meynetetal. 1994 ) as observed in 7^ Vel. Therefore, these 
models would predict a minimum -^Al yield of 5.5 x 10^^ M0 
for 7^ Vel with initial mass of the WR component Ali = 40 M0 
and even 1.2 x 10"* M© for Mi = 60 Mq, the value used by 



Meynetet al. (1997 ) for the description of 7^ Vel. These yields 
are co nsistent with the yiel d s predicted by the sing le-star mod- 
els of Langer et al. (1995 ). Schaerer et al. (1997 ) estimate an 
initial mass of Mi = 57 ± 15 M0 for the WR star, using the 
Geneva models. 

Has all ^^Al been expelled yet or is some fraction still 
buried invisibly inside the star? The spectral type of the Wolf- 
Rayet star is WC, i.e. the stellar wind is carbon rich, which 
means that products of core helium burning have reached the 
stellar surface. Since remaining ^^Al in the stellar core is ef- 
ficiently destroyed in helium burning due to neutron captures, 
standard stellar evolution models predict that the wind ejection 
phase of hydrostatically produced '^ Al has ended already. 

How much ^*A1 has already decayed? Typical WR lifetimes 
from observational constraints and theoretical models are on 
the order of Si 5 x 10^ yr, which is shorter than the ^^Al lifetime 
by a factor of 2. This could account for a reduction of observ- 
able ^''Al in the ISM of at most 40% if the WR star were close 
to the end of its evolution. Some models predict WR lifetimes 
in excess of a million years for the most massive stars (Meynet 
et al. 1997) but these models also predict sufficiently large ^^Al 
yields such that the prolonged time available for decay would 
not reduce the observable amount below the measured limit. 
Also note that WR lifetimes strongly depend on the mass loss 
description applied in the model. 

A straight-forward explanation for the missing 1.8 MeV 
flux would be sub-solar metallicity since the ^^Al yields scale 
approximately like Z^. Yet, sub-solar abundances in the ISM 
observed in the Vela direction can readily be understood by 
dust formation (and therefore depletion of the gaseous phase) 



rather than by intrinsic low metallicity (Fitzpatrick 1996) 
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B: p=0 (BL95,L95) a 
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Fig. 3. ^^ Al yields as a function of the initial mass of the WR 
component in the binary system. Plotted are single star ("S") 
and binary ("B") models as well as models of rotating (sin- 
gle) stars ("R"), all at solar metallicity. L95: Langer et al. 



(1995); M97: Meynet et al. (1997); BL95: Braun & Langer 



(1995); L97: Langer et al. (1997). For interpolation, the sin 



Analysis of spectroscopic data for the O star in the binary sys- 
tem is indeed consistent with solar metallicity (De Marco & 
Schmutz 1999). 



gle star models were fitted with a power law (dashed and dash- 
dotted lines). The horizontal solid line marks our 2 a limit (us- 
ing model c) for the ^^Al mass surrounding the WR star The 
light-gray region corresponds to the 1 a distance uncertainty, 
the dark-gray region is excluded by the 1.8 MeV flux limit. 
The high-lighted area indicates the most probable initial mass 
range. Single star models are barely compatible with our mea- 
surement, and even binary models are consistent only at the 
lowest allowed initial mass range and would require significant 
mass transfer (/3) from the primary to the secondary. 

With these considerations. Fig. shows that ^^Al yields 
predicted from single star models are barely compatible with 
the COMPTEL flux limit. This would suggest that model pa- 
rameters such as the mass loss description, which has great- 
est impact on the stellar structure for initial masses ^ 40 Mq, 
or internal mixing parameters like core overshooting or semi- 
convection may have to be modified. 

3.4. Binary models 

Differences in the stellar evolution of the primary star in a rela- 
tively wide binary system such as 7^ Vel (a « 1 AU) stem from 
mass loss, while effects of tidal forces on the stellar structure 
are negligible during the main sequence phase when ^''Al is 
produced in the core. In addition to the mass loss mechanisms 
of single stars, Roche Lobe Overflow (RLOF) can change the 
stellar structure of both primary and secondary star For the 
discussion of ^^Al yields, we can concentrate on the evolu- 
tion of the primary since the secondary O star merely reveals 
unprocessed material from the stellar envelope at its surface, 
but might bur y some of the p rocessed material transfered from 
the primary. Langer (1995 ) proposes that binary stars with 
Mi ^ 40 Mq undergoing RLOF are essentially indistinguish- 



U. Oberlack et al.: COMPTEL Limits on -*A1 1.809 MeV line emission from 7^ Vel 





8 


9 


10 11 


12 


13 


1 


1.0 












' ' _ 


0.8 


r 










-_ 


0.6 


r ^ 


^.^ 








-_ 


0.4 


- 










-_ 


0.2 


- 










— :; 


0.0 












'- 



40 Mq the lowest possible initial mass for the WR star, the 
fact that binaries with initial mass ratios (secondary / primary) 
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Fig. 4. Maximum value for the mass fraction (3 transferee! to 
the secondary assuming current masses Mwr — 9 Mq, Mq — 
29 Mq, and an initial m ass ratio q > 0.2 to avoid merging 
dVanbeveren et al. 19981). 



able from single stars in the same mass range undergoing a 
phase of very intense mass loss as LBVs. Differences in the 
^^Al yield in the circumstellar medium would merely result 
from the fraction of expelled mass captured by the secondary 
companion, i.e. some fraction S5 (3 could be buried in the sur- 
face layers of the secondary. (Stellar winds from primary and 
secondary will always transport some fraction of ^^Al into the 
ISM.) Vanbeveren (1991) argues that this LBV mass loss may 
even prevent any occurence of RLOF, which means that ^*A1 
yields should remain the same as for single stars. 

For initial masses below ^ 40 Mq, RLOF can provide 
additional mass loss not attainable in the single star scenario, 
pushing down the lower i nitial mass limit for the formation of 
WR stars to about 20 Mq ( |Vanbeveren et al. 1998| ). Yet, adopt- 
ing current masses of 9 Mq for the WR star and 29 Mq for the 
O star and assuming a fraction of mass transfered to the sec- 
ondary as high as (3 = 0.5 would yield a minimum initial mass 
for the primary of ~ 30 Mq, given that the primary star must 
have been the more massive partner initially to evolve faster 
The observed current mass loss rate of (2.8^0 g) x 10^^ MQ/yr 
( Schaerer et al. 1997 ), quite typical for this type of star, sup- 
ports a minimum mass lost into the ISM of at least 10 Mq 



in the last few 100,000 years. Vanbeveren et al. (1998 ) quote 
a minimum initial mass of 38 Mq for the WC star in 7^^ Vel 
based on detailed models of binary evolution. 

Overall, the minimum initial mass for the WR star in bi- 
nary models is found close to the values obtained by single 
star models, namely around 40 Mq. While the real initial mass 
may well have been larger, it is apparent from Fig. Wthat dis- 
crepancies between predicted ^^Al yields and the measured 
1.8 MeV flux limit become quite severe for larger masses. Even 
for il/j = 40 Mq, models are in clearly better agreement with 
the flux limit if a significant fraction of the ejected mass of 
the primary accreted onto the secondary. The binary models 
of Braun & Langer (1995) lead to a typical reduction in ^''Al 
yield of about 40% for models with (3 — 0.5 as displayed in the 
figure. 

Information regarding a quantitative estimate of (3 is still 
sparse, but some statements can be made. If we consider il/j « 



q < 0.2 are expected to merge (Vanbeveren et al. 1998) leads 
to a minimum initial mass for the O star of ~ 8 Mq, hence 
an initial total mass of the system of <; 48 Mq. Assuming a 
current mass of (29 + 9 = 38) Mq as for the initial mass es- 
timate, a minimum of 10 Mq must have been expelled to the 
ISM, while about 20 Mq would have been transfered to the 
secondary, corresponding to /5 w 2/3. Any larger value of q, 
i.e. larger initial mass of the O star, would yield lower (3. For 
larger values of the initial mass of the WR star, the maximum 
allowed /3 for q > 0.2 rapidly decreases as illustrated in Fig. ^ 
Observation of an "ejecta-type" ring nebula around 7^ Vel 



with 57' angular diameter (Marston et al. 1994), correspond 



ing to a radius of 2.1 pc, demonstrates that significant mass has 
been expelled from the system, probably during a preceding 
LBV / WN phase, even though the mass of the shell has not yet 
been determined. If indeed the IRAS Vela shell were the rem- 
nant of the main sequence bubble of the primary star, the mass 
transfer to the secondary would be negligible compared to the 
ejected mass and the WR star could be reasonab ly modeled by 
a Mi « 60 Mq star ( Garcia-Seguraet al. 1996 ). This scenario 
is supported by the recent finding that the helium abundance in 
the companion O star is not enriched (De Marco & Schmutz 
1999) as one might expect if significant amounts of processed 
material had been transfered to the secondary star. This, how- 
ever, would imply a major conflict of the ^^Al yields predicted 
by corresponding stellar models with our measurement. 

There are additional qualitative arguments for both very lit- 
tle mass transfer to the secondary but also some mass trans- 
fer: The relatively high excentricity of the orbit favours little 
mass transfer, which usually tends to reduce the excentricity. 
On the other hand, the high rotation velocity of the secondary 
ofurotsini — 220 kms^^ (Baadeetal. 1990) may be the result 
of spinning up by accretion. More detailed modeling is required 
to resolve these issues. 



4. Conclusion 

Given the small distance of 7^ Vel from HIPPARCOS measure- 
ments and predicted ^^Al yields of current stellar models, the 
non-detection of 1.8 MeV emission by COMPTEL comes as a 
surprise. Combined with other observations regarding current 
masses, metallicity and mass transfer, only a very small vol- 
ume of the model parameter space remains consistent with the 
COMPTEL 2 crupper limit of A/WRii < (6.3^1:4) IQ-^Mq. 
Single star models are in conflict with this value. Binary mod- 
els alleviate the discrepancy if significant mass transfer to the 
secondary occured, burying some fraction of ^^Al in the surface 
layers of the O star, and if the initial mass of the WR star was 
close to its minimum value of ~ 40 Mq. It may be more likely, 
however, that adjustment of some of the model parameters, e.g. 
the parametrization of mass loss, is required. 

Unfortunately, 7^ Vel is the only known WR star for which 
there was hope to obtain a positive detection with current in- 
struments. The WR stars next closest to the Sun (WR142, 
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WR145, and WR147 in the Cygnus region) are at least a fac- 
tor of two more distant, and therefore their expected fluxes 
are out of reach for COMPTEL. For the forthcoming INTE- 
GRAL mission, these stars may just become detectable, and the 
1.8 MeV flux from 7^ Vel will be tested down to significantly 
lower values. Yet, to probe the 1.8 MeV flux from individual 
WR stars within the radius of completeness of the current WR 
catalogue (~ 2.5-3 kpc), would require a next-generation in- 
strument with a sensitivity of 10~^ 7 cm^^ s^^ and an angular 
resolution < 0.2°. 
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